Maximally entangled mixed states: Creation and concentration 
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Using correlated photons from parametric downconversion, we extend the boundaries of experi- 
mentally accessible two-qubit Hilbert space. Specifically, we have created and characterized maxi- 
mally entangled mixed states (MEMS) that lie above the Werner boundary in the linear entropy- 
tangle plane. In addition, we demonstrate that such states can be efficiently concentrated, simul- 
taneously increasing both the purity and the degree of entanglement. We investigate a previously 
unsuspected sensitivity imbalance in common state measures, i.e., the tangle, linear entropy, and 
fidelity. 
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By exploiting quantum mechanics it is possible to 
implement provably secure cryptography jlj, teleporta- 
tion 13] , and super-dense coding . These protocols and 
most others in quantum information processing require a 
known initial quantum state, and typically have optimal 
results for pure, maximally entangled initial states. How- 
ever, decoherence and dissipation may cause the states 
to become mixed and/or less entangled. As the suc- 
cess of a protocol such as quantum teleportation often 
hinges on both the purity and the entanglement of the 
initial state 0, it is important to study the interplay of 
these properties. Using a source of 2-qubit polarization 
states , we investigate the creation of maximally en- 
tangled mixed states, and their concentration [liQjSH- 

Entangled states have been demonstrated in a variety 
of systems [H Ell E H H [13. In fact, there are 
several classes of entangled states; maximally entangled 
and nonmaximally entangled pure states [3, lli ll^ > non- 
maximally entangled mixed states [T^ . and the special 
case of Werner states ^3 (incoherent combinations of a 
completely mixed state and a maximally entangled pure 
state) have all been experimentally realized using optical 
qubits. For some time it was believed that Werner states 
possess the most entanglement for a given level of mixed- 
ness. However, Munro et al. ^3 discovered a class of 
states that are more entangled than Werner states of the 
same purity. These maximally entangled mixed states 
(MEMS) possess the maximal amount of entanglement 
(tangle or entanglement of formation) for a given degree 
of mixedness (linear entropy) [20II2H . 

By generating states close to the MEMS boundary, 
we have experimentally explored the region above the 
Werner state line on the linear entropy-tangle plane [23 ■ 
We have also implemented a partial-polarizer filtra- 
tion/concentration technique which simultaneously in- 
creases both purity and entanglement, at the cost of de- 
creasing the ensemble size of initial photon pairs. Though 
the implementation requires initial state knowledge, we 
show that MEMS exist for which this "Procrustean" fil- 
tering technique 0, H, E3 is much more efficient than 
other recent entanglement concentration schemes [l^, 



|25|, even after modification to work on MEMS. 

The exact form of the MEMS density matrix de- 
pends on the measures used to quantify the entan- 
glement and mixedness [2li |: here we use the tangle 
(r(p) = [max{Q, Al -A2- A3 -A4}]^) 01, i-e., the con- 
currence squared; and the linear entropy [Sl {p) = f [1 ^ 
Tr(p^)]) y. Here Ai are the square roots of the eigen- 
values of p{j72 ® (T2)p*{a'2 (8) CT2), in non- increasing order 

by magnitude, with (T2 = ^ ^ 0^ ) ' ^^^^ parame- 
terization, where r is the concurrence, the MEMS den- 
sity matrices exist in two subclasses yJl, Pmems i and 
Pmems //, which have two and three eigenvalues, respec- 
tively: 
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Our creation of MEMS involves three steps: creating 
an initial state of arbitrary entanglement, applying lo- 
cal unitary transformations, and inducing decoherence. 
First, frequency degenerate 702-nm photons are created 
by pumping thin nonlinear /3-Barium Borate (BBO) crys- 
tals with a 351-nm Ar-ion laser. Polarization entangle- 
ment is realized by pumping two such crystals oriented 
such that their optic axes are in perpendicular planes. 
With a pump polarized at 9i, a variable entanglement 
superposition state cos6i\HH) + sm9i\VV) is created, 
where \HH) represents two horizonta-lly polarized and 
\VV) two vertically polarized photons [5|,|l6j- The pump 
polarization is controlled by a half- wave plate (HWPi in 
Fig. HI set to 6*1/2. 

To create the MEMS I, we start by setting the ini- 



2 



HWPi BBO 



SPDC 
Source 




MEMS 
^Preparation^ 



State 
iConcentration; 



HWP IRIS LENS 

Tomography System 



FIG. 1: Experimental arrangement to create, and concen- 
trate MEMS. A half-waveplate (HWPi) sets the initial en- 
tanglement of the pure state. The i^-plate sets the relative 
phase between \HH) and \VV) in the initial state. HWP2 
and HWP3 rotate the state into the active bases of the de- 
coherers to adjust the amount of entropy. The tomography 
system uses a quarter- waveplate (QWP), HWP, and a polar- 
izer in each arm to analyze in arbitrary polarization bases; the 
transmitted photons are counted in coincidence via avalanche 
photodiodes. The dashed box contains HWP4 (oriented to 
rotate \H) ^ \V) in the first arm of the experiment) and 
concentrating elements (a variable number of glass pieces ori- 
ented at Brewster's angle to completely transmit \H), but 
only partially transmit |V^)). 

tial degree of entanglement to that of the tar get MEMS. 
Next a maximum likelihood tomography [lil l27l | of this 
initial entangled state is taken and used to numerically 
determine the appropriate settings of HWP2 and HWP3 
in Fig. n These waveplates set the diagonal elements of 
the density matrix to the target values for the desired 
MEMS. The initial tomography must be precise, because 
the waveplate settings are critically dependent on the ini- 
tial state, as well as on the precise birefringent retarda- 
tion of the waveplates themselves. After the waveplates, 
the state passes through decoherers, which lower specific 
off-diagonal elements in the density matrix, yielding the 
final state. In our scheme, each decoherer is a thick bire- 
fringent element (^1 cm quartz, with optic axis horizon- 
tal) chosen to have a polarization-dependent optical path 
length difference ('--^140A greater than the downcon- 
verted photons' coherence length (Lc = A^/AA ^ 70A, 
determined by a 10-nm FWHM interference filter placed 
before each detector), but much less than the coherence 
length of the pump [23 . 

The decoherer in each arm couples the polarization 
with the relative arrival times of the photons [s^l . While 
two horizontal {\HH)) or two vertical (|1^^)) photons 
will be detected at the same time, the state \HV) will 
in principle be detected first in arm one and then in arm 
two, and vice versa for the state l^-ff) (assuming the 
decoherer slows vertically polarized photons relative to 
horizontally polarized ones). Tracing over timing infor- 
mation during state analysis then erases coherence be- 
tween any distinguishable terms of the state (i.e., only 
the coherence term between \HH) and \VV) remains). 
A sample tomography of a MEMS I is shown in Fig.l^^a). 

MEMS II are created by first producing the MEMS I 
at the MEMS I/II boundary, i.e., the state with = |. 
In order to travel along the MEMS II curve, the optical 
path length difference in one arm must be varied from 
140A. This couples different relative timings to the \HH) 
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FIG. 2: MEMS data, (a) Density matrix plots of the real com- 
ponents of a MEMS I (r = I) and a MEMS II (r = 0.3651). 
The imaginary components are negligible (on average less 
than 0.02) and not shown, (b) Linear entropy-tangle plane. 
Shown are theoretical curves for MEMS I (solid hue) , MEMS 
II (dashed line), and Werner states (dotted line). Four target 
MEMS are indicated by stars; experimental realizations are 
shown as squares with error bars. The shaded patches around 
each target state show the tangle (T) and linear entropy {Sl) 
for 5000 numerically generated density matrices that have at 
least 0.99 fidelity [sj with the target state. T^O (1) corre- 
sponds to a product (maximally entangled) state. Sl=G (1) 
corresponds to a pure (completely mixed) state. 

and \VV) states, reducing the coherence between them. 
For instance, to make the MEMS II (B) in Fig. \^&), 
140A decoherence was used in one arm, 90A in the other. 
Fig. IHa) indicates very good agreement between the ory 
and experiment with fidelities of ~99% (the fidelity |31| 
between the target state pt and the measured state pm is 
given by F{pt,pm) = |Tr {.J ^tPrn^/pt) \ )■ 

The states (A) and (B) are shown in the Sl-T plane 
in Fig. IHb), along with other MEMS we created. The 
states do not hit their Sl-T targets (shown as stars in 
the figure) within errors, even though the states have 
very high fidelities (> 99%) with their respective tar- 
gets. To explore the discrepancy, for each target we nu- 
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FIG. 3: Concentration data. Shown are concentrations for 
three initial states, A (triangles) and B (filled squares) as in 
Fig. 121 and C (open squares), along with the number of partial 
polarizing glass pieces in each arm. The expected concen- 
trated state path, calculated using 01, is shown with stars. 
The concentrated states agree with theory for small numbers 
of glass pieces, but as more slips are used, the state concen- 
trates better than expected. We believe this is due to the 
extreme sensitivity of the trajectory to small changes in the 
initial state. However, even in theory, excessive filtration will 
eventually produce a pure product state (shown as an exten- 
sion of A's theory curve), due to small errors in the initial 
MEMS. 

merically generated 5000 density matrices that had at 
least 0.99 fidelity with the target density matrix. The 
Sl and T of the numerically generated states are plotted 
in Fig. E^b) as shaded regions surrounding the targets. 
The fact that these regions are rather large (and overlap 
with our measured MEMS) explains our results, but is 
surprising nonetheless. The unexpectedly large size of 
these patches arises from the great difference in sensi- 
tivity between the state measures of fidelity, tangle and 
entropy: for small perturbations {5r) of the MEMS pa- 
rameter r, the fidelity is only quadratic in 6r, while Sl 
and T are linear in Sr [3^ . 

While our initial goal was to produce states of maxi- 
mal tangle for a given linear entropy, maximally entan- 
gled pure states are generally more useful for quantum 
information protocols. However, in some cases, weakly 
entangled mixed states may be the only available re- 
source. It is therefore important to study ways to si- 
multaneously decrease the entropy and increase the en- 
tanglement of an ensemble of photon pairs (necessarily 
at the cost of reducing the size of the ensemble). Re- 
cently several such entanglement concentration experi- 
ments have been reported, relying on two-photon inter- 
ference effects ll^, . An interesting characteristic of 
MEMS is that they can be readily concentrated by a 
"Procrustean" method of local filtering . To con- 

centrate we first modify the MEMS using IIWP4 at 45° to 
exchange \H) <-!■ \V) in the first arm, changing the non- 
zero diagonal elements of the MEMS density matrix to 
\HV){HV\, \VH){VHl and \VV){VV\. By reducing the 
element of the rotated MEMS, the outcome 



will be driven toward the maximally entangled pure state 
|(/)+) = i\HV) + \VH))/V2. We achieve this by insert- 
ing glass pieces (each piece consisting of four '^Imm thick 
microscope slides sandwiched together with index match- 
ing fiuid) oriented at Brewster's angle, as indicated in the 
dotted box in Fig.^ Equal numbers of pieces are used in 
both arms; they are oriented to nearly perfectly transmit 
horizontally polarized photons (transmission probability 
Th = 0.990 ± 0.006) while partially reflecting vertically 
polarized photons (Ty = 0.740 ± 0.002). 

We concentrated a variety of MEMS. Fig. shows the 
results for the MEMS I and II of Fig. Inland an additional 
MEMS I (C). As the number of glass pieces is increased, 
the states initially become more like a pure maximally en- 
tangled state. For example, in the case of (A), the fidelity 
of the initial MEMS with the state |(/)+) is 0.672. When 
the state is concentrated with eight glass slips per arm, 
the fidelity with is 0.902; 4.5% of the initial photon 
pairs survive this filtering process. The theoretical max- 
imum survival probability is 6.4%. Note a characteristic 
difference between the two MEMS subclasses: MEMS II 
cannot be filtered into a Bell state. 

We now compare the theoretical efficiency of our lo- 
cal filtering scheme with the interference-based concen- 
tration proposal of Bennett et al. 0, assuming identical 
initial MEMS and the same number of photon pairs. We 
shall compare the average final entanglement of forma- 
tion (Ep) 2^ (i.e., the Ep of the concentrated state mul- 
tiplied by the probability of success) per initial pair. The 
Bennett el al. 0] scheme was recently approximated by 
Pan et al. [i^l, with CNOT operations replaced by po- 
larizing beam splitters; however, due to incomplete Bell 
state analysis, the probability of successful concentration 
is only 50% of the original proposal (the recent scheme 
of Yamamoto et al. [23 is unable to distill MEMS). The 
first ste p o f both schemes is to perform a "twirling" op- 
eration [33 to transform a general entangled state into a 
Werner state. However, this initial operation usually de- 
creases the entanglement, and the scheme with twirling 
is efficient only when r is close to 1. 

In fact, MEMS I could also be distilled without the 
twirling operation, using the scheme of Pan et al. But 
then the probability of success depends on the parame- 
ter r. For most MEMS, the maximum distillation effi- 
ciency from filtration can exceed that achievable using 
the interference-based methods [s^. For example, as 
shown in Table I, when the initial state is a MEMS with 
r = 0.778, the two-piece filtering technique has a theo- 
retical Ep per pair nearly three times higher than the 
interference scheme without twirling, even though a suc- 
cessful concentration produces nearly the same Ep. In 
theory, using 2 to 5 slips achieves both higher entangle- 
ment of the successful state and better average entangle- 
ment yield. In practice, the filtration technique is much 
more efficient (see the final columns of Table I) [s^ . 

We have demonstrated a tunable source of high fidelity 
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14.2% 


0.93 


0.13 
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TABLE I: EfRciency comparison of concentration technique 
of Bennett et al. using ideal CNOT 0, interference-based 
concentration [2^ without twirling, and Procrustean filter- 
ing, for an initial MEMS with r = 0.778 and Ef = 0.69. The 
scheme of Bennett et al. requires a twirling operation that de- 
creases the initial Ep to 0.418 before the concentration f3^. 
In all schemes, except for the final column, we assume the 
ideal case, i.e., no loss and perfect detector efficiency. To cal- 
culate the no-loss result for our filtering scheme, we normalize 
the measured partial polarizer transmission coefficients (of a 
single glass piece) to Th = 0.740/0.990 and Ty = 1. In the 
interference schemes, columns 2-4 assume the existence of the 
required two identical pairs, but in practice this requirement 
is difficult to achieve |34| . This limitation is refiected in col- 
umn 5, which lists the average Ef per initial pair achieved in 
our experiment, to be compared with the much lower value 
achievable with current interference method technology. 

MEMS. As a consequence of comparing the T-Sl and 
fidelity values of generated MEMS with the theoretical 
targets, we identify and explain an unsuspected difference 
in sensitivity in these state measures. Furthermore, we 
have applied a Procrustean filtering technique to several 
MEMS, realizing a measured efficiency that is well above 
that achievable using other methods. However, in the 
limit of very strong filtering, small perturbations in the 
initial state will eventually dominate the process, yielding 
only product states (see Fig. 3). In practice, therefore, 
it may be optimal to combine both methods. 
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1 + ^5r, and f (p,(r„), p,(r)) 



r(p.(r))/T(p,(ro)) 
1 + Bi{Sr)^, where the subscript i denotes the class of 
MEMS, and the constants are given by Ai — ^r"""*" , , 
An = Bi = -1 ™^ f ■- - 3 



4ro(l-ro) 



, and Bii = 



2(9r-i^-4) • 



In "twirling," a random SU(2) rotation is independently 
performed on each photon pair. 
[34] Because it is presently very difficult to produce simul- 
taneous indistinguishable pairs of photons, the filtration 
technique is much more efficient, e.g., where typically 
20% of our incident ensemble of pairs survived, less than 
0.005% would survive (estimated from the 2-fold and 4- 
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fold coincidence data reported in |2^) in the interference 
schemes, which require 4 photons. 



